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ABSTRACT 

We measure the evolution of the luminous red galaxy (LRG) luminosity function in the redshift range 
0.1 < z < 0.9 using samples of galaxies from the Sloan Digital Sky Survey as well as new spectroscopy 
of high-redshift massive red galaxies. Our high-redshift sample of galaxies is largest spectroscopic 
sample of massive red galaxies at z ~ 0.9 collected to date and covers 7 deg^, minimizing the impact 
of large scale structure on our results. We find that the LRG population has evolved little beyond the 
passive fading of its stellar populations since z ~ 0.9. Based on our luminosity function measurements 
and assuming a non-evolving Salpeter stellar initial mass function, we find that the most massive 
(L > 3L*) red galaxies have grown by less than 50% (at 99% confidence), since z = 0.9, in stark 
contrast to the factor of 2-4 growth observed in the L* red galaxy population over the same epoch. 
We also investigate the evolution of the average LRG spectrum since z 0.9 and find the high-redshift 
composite to be well-described as a passively evolving example of the composite galaxy observed at 
low-redshift. From spectral fits to the composite spectra, we find at most 5% of the stellar mass in 
massive red galaxies may have formed within IGyr of z = 0.9. While L* red galaxies are clearly 
assembled at z < 1, 3L* galaxies appear to be largely in place and evolve little beyond the passive 
evolution of their stellar populations over the last half of cosmic history. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: photometry - 
galaxies: statistics - galaxies: fundamental parameters 



1. INTRODUCTION 

The favored model for the evolution of galaxies is 
through the hierarchical merging of smaller satellite 
galaxies into larger systems. The details of the frequency 
and efficiency of the merging process are poorly con- 
strained, especially in the densest environments. As the 
endpoint of the hierarchical merging process, the most 
massive galaxies are most sensitive to various merger 
models assumptions and thus offer a strong opportunity 
to constrain models of galaxy formation and evolution. 

Observations of the evolution of early-type galaxy stel- 
lar populations have shown that the stars in these galax- 
ies formed at z > 2 and that the galaxies have had little 
star formation since that epoch (Bower et al. 1992; Ellis 
et al. 1997; Kodama et al. 1998; de Propris et al. 1999; 
Brough et al. 2002; Holden et al. 2005; Wake et al. 2005; 
Pimbblet et al. 2006; Jimenez et al. 2006; Bernardi et 
al. 2003a,b,c,d; Glazebrook et al. 2004; McCarthy et al. 
2004; Papovich et al. 2005; Thomas et al. 2005; Bernardi 
et al. 2006; Cool et al. 2006). While the average popula- 
tion of massive galaxies appears to be quite old and pas- 
sively evolving, a number of studies have indicated that 
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local massive early-type galaxies show signs of recent star 
formation activity (Trager et al. 2000; Goto et al. 2003; 
Fukugita et al. 2004; Balogh et al. 2005). The fraction 
of early-type galaxies with evidence of recent star forma- 
tion seems to increase to high redshift and decreases with 
increasing stellar mass (Le Borgne et al. 2006; Caldwell 
et al. 2003; Nelan et al. 2005; Clemens et al. 2006). 

At z < 1, early- type galaxies form a tight relationship 
between their rest-frame color and luminosity (the so- 
called color-magnitude relation or red-sequence of galax- 
ies) wherein more luminous (and hence more massive) 
galaxies have redder colors then less-massive counter- 
parts (Visvanathan & Sandage 1977; Bower et al. 1992; 
Hogg et al. 2004; Mcintosh et al. 2005; Willmer et al. 

2006) . The tight dispersion around this relationship 
implies that, at fixed luminosity, galaxies on the red- 
sequence share very similar star formations histories. If 
massive galaxies have undergone any mergers since z ~ 1 , 
the mergers must have resulted in very little star forma- 
tion; the addition of even a small fraction of blue stars 
would result in a larger intrinsic scatter than observed 
(Cool et al. 2006). 

The extent to which gas-poor mergers that result in no 
new star formation are involved in the build-up of mas- 
sive galaxies is a topic of much current research. While 
examples of these mergers have been observed at low 
redshift (Lauer 1988; van Dokkum 2005; Mcintosh et al. 

2007) and at intermediate redshifts (van Dokkum et al. 
1999; Bell et al. 2006b; Tran et al. 2005; Rines et al. 2007; 
Lotz et al. 2008), the extent to which massive galaxies 
participate in these merger events is controversial. Bell 
et al. (2006a) and Le Fevre et al. (2000) estimate that L* 
red galaxies experience 0.5-2 major mergers since z ~ 1.0 
based on pair counts of galaxies, van Dokkum (2005) 
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identified galaxies which have hkely undergone a recent 
gas-poor merger based on the presence of diffuse emis- 
sion extended from the main galaxies and found that 
35% of today's bulge dominated galaxies have experi- 
enced a merger with mass ratio greater than 1:4 since 
z ^ 1. Based on the very small-scale correlation func- 
tion of luminous red galaxies from SDSS, Masjedi et al. 
(2006) concluded that mergers between these very mas- 
sive systems occur quite rarely at z ~ 0.3 with rates 
< 1/160 Gyr"\ Masjedi et al. (2007) calculate that mas- 
sive early- type galaxies have grown by 1.7% per Gyr on 
average since z ~ 0.2 due to mergers with all other galax- 
ies. 

Studies based on the number counts of galaxies from 
COMBO-17, DEEP2, and the NOAO Deep Wide-Field 

Survey (NDWFS) all agree that the stellar mass aver- 
aged all red-galaxies has at least doubled since z ~ 1 
(Brown et al. 2007; Willmer et al. 2006; BeU et al. 2004). 
While the truncation of star formation in blue galaxies 
and subsequent passive fading of the stellar populations 
can explain the growth of L* galaxies since z ^ 1, the 
lack of very massive blue galaxies at redshift of unity 
(Bell et al. 2004) indicates that any evolution of the most 
massive galaxies must be fueled by mergers of less lumi- 
nous red-galaxies and not from pure passive evolution of 
massive star forming galaxies. While red galaxies with 
L K, L* appear to grow substantially since z ~ 1, re- 
sults from Brown et al. (2007) indicate that very lumi- 
nous {L > 4L*) galaxies have grown by only 25% since 
z ~ 1.0. Similarly, Wake et al. (2006), used a combina- 
tion of the SDSS and 2dF-SDSS LRG and QSO (2SLAQ) 
sample to measure the evolution of the massive galaxy 
luminosity functions to z = 0.6 and found that at least 
half of the massive early-type galaxies present at z = 0.2 
must have been well assembled by z ~ 0.6. These inves- 
tigations agree with a number of studies which have sug- 
gested little or no evolution in the most massive galaxy 
populations (Lilly et al. 1995; Lin et al. 1999; Chen et al. 
2003; Bundy et al. 2006; Willmer et al. 2006; Cimatti et 
al. 2006). 

In this paper, we present new observations of massive 
red galaxies at 0.7 < z < 0.9 and augment it with sam- 
ples of massive red-sequence galaxies from SDSS in order 
to quantify the evolution of the massive galaxy luminos- 
ity function over half of cosmic history. Our high-redshift 
spectroscopic survey is unaffected by possible systematic 
errors from photometric redshifts and covers 7 square 
degrees, minimizing the effects of cosmic variance due to 
large-scale galaxy clustering. 

After describing our galaxy sample selection criteria in 
§2, we discuss the construction of our massive red galaxy 
luminosity functions in §3. In §4, we interpret out lumi- 
nosity function measurements and examine the compos- 
ite spectrum of massive red galaxies since z ~ 0.9 in §5 
before closing in §6. All magnitudes discussed in the text 
are AB (Oke & Gunn 1983). When calculating luminosi- 
ties and volumes, wc use the cosmological world model of 
0„ = 0.25, n^n + ^A = 1, and Ho = 100 h km s^^ MpcT^ 
(Spergel et al. 2007). When calculating time, for exam- 
ple when considering the aging of stellar populations, we 
use h = 0.7. All magnitudes are corrected for dust ex- 
tinction using the dust maps of Schlegel, Finkbeiner, & 
Davis (1998). 



2. SAMPLE CONSTRUCTION 

2.1. SDSS Galaxy Sample 

The Sloan Digital Sky Survey (SDSS; York et al. 2000; 
Adelman-McCarthy et al. 2007) has imaged tt stera- 

dians of the sky in five bands, ugriz, (Fukugita et al. 
1996) with a dedicated 2.5m telescope located at Apache 
Point Observatory (Gunn et al. 2006). Imaging is per- 
formed with a CCD mosaic in drift-scan mode (Gunn 
et al. 1998) with an effective exposure time of 54s. Af- 
ter images are reduced (Lupton et al. 2001; Stoughton 
et al. 2002; Pier et al. 2003) and calibrated (Hogg, et 
al. 2001; Smith et al. 2002; Ivezic et al. 2004; Tucker et 
al. 2006), objects are chosen for follow-up spectroscopy 
using an automated spectroscopic fiber assignment algo- 
rithm (Blanton et al. 2003a). Two galaxy samples are se- 
lected for spectroscopy from SDSS imaging. The MAIN 
galaxy sample (Strauss et al. 2002) is a complete, flux- 
limited (r < 17.77), sample of galaxies with an average 
redshift of 0.1. The Luminous Red Galaxy (LRG) sam- 
ple (Eisenstein et al. 2001) selects luminous early-type 
galaxies out to z ~ 0.5 with r < 19.5 using several color- 
magnitude cuts in g, r, and i. The average redshift of 
the LRG sample is ^ 0.3. 

In addition to its contiguous coverage of the Northern 
Galactic cap, the SDSS also conducts a deep imaging 
survey, SDSS Southern Survey, by repeatedly imaging 
an area on the celestial equator in the Southern Galactic 
Cap. The data we utilize here includes 300 deg^ of imag- 
ing that has been observed an average of 20 times and 
up to 30 times. Objects detected in each observational 
epoch were matched using a tolerance of 0.5 arcseconds 
to create the final coadded catalog. The measured pho- 
tometry from each epoch were combined by converting 
the reported asinh magnitudes (Lupton et al. 1999) to 
flux and then calculating the mean value. Errors on each 
parameter are reported as the standard deviation of the 
flux measurements. 

While the LRG color selection criteria identify mas- 
sive red galaxies at moderate redshifts, at redshifts below 
z ~ 0.2 the LRG color selection becomes too permissive - 
under-luminous blue galaxies are allowed into the sample 
(Eisenstein et al. 2001). In order to construct a sample 
of galaxies at 0.1 < z < 0.2, we thus rely on the MAIN 
galaxy sample; in this redshift range, the massive galax- 
ies of interest pass the r < 17.77 flux hmit of the MAIN 
sample. We utilize a simple rest-frame color- luminosity 
cut, Mg < — 21 and {g — i)rost > 2 to select low-redshift 
galaxies on the red-sequence. These cuts result in 23,854 
LRGs at O.K z < 0.2. At 0.2 < z < 0.4, the LRG selec- 
tion provides a clean sample of 46,856 massive red galax- 
ies which we consider our intermediate redshift galaxy 
sample. Our low- and intermediate- redshift samples 
clearly have quite different selection functions in their 
rest-frame colors which immt be considered when mea- 
suring the evolution between samples; we will address 
this when we present our luminosity function measure- 
ments in §3. 

2.2. SDSS Photometry 

As described in detail in Stoughton et al. (2002), 
Strauss et al. (2002), and Blanton et al. (2001), SDSS 
galaxy photometry is reported using two systems. Each 
galaxy in SDSS is flt by two seeing-convolved models, 
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a pure do Vaucouleurs (1948) model and a pure expo- 
nential profile. The best-fitting model in the r-band is 
used to determine the flux of the galaxy in each of the 
other bands by adjusting the normalization to the model 
while leaving all other parameters fixed to those derived 
in the r-band. Alternatively, the Petrosian magnitude is 
defined to be the flux within 29 p where 9p is defined to 
be the radius at which point 

_ J^:^^S d9'2^e'm I [7r(1.252 - 0.852)^2] 
\l dr'2TT9'I{9') / [7r6i2] 

(1) 

falls below 0.2. Here, I{9) is the azimuthally averaged 
surface brightness profile of the galaxy. The Petrosian ra- 
dius is determined in the r-band and then applied to each 
of the other bands. While the Petrosian flux measure- 
ment contains a constant fraction of the galaxy's light 
in the absence of seeing, independent of its size or dis- 
tance, model magnitudes are unbiased in the absence of 
color gradients and provide a higher signal-to-noise ratio 
color measurement than Petrosian colors. As the Pet- 
rosian flux aperture is defined based on the shape of 
the light distribution, it doesn't require measuring the 
faint, low-surface brightness, isophotes of the galaxy at 
large radius which is quite difficult with shallow photom- 
etry. Throughout this paper, we use model magnitudes 
when discussing colors of galaxies and Petrosian quanti- 
ties when calculating luminosities. 

As has been noted by Lauer et al. (2007), SDSS pho- 
tometry of very large (rdr > 10") galaxies at low red- 
shift have large systematic differences from measured 
photometry in the literature. For very large galaxies, 
the automated photometric pipeline includes galaxy light 
in the estimation of the local sky background and thus 
underestimates the total galaxy flux. At z > 0.1, we 
expect this effect to play a minimal role and thus per- 
form no correction to our photometry. In order to en- 
sure that this is a valid approach, we simulate 2,000 
galaxies at 0.1 < z < 0.4 with properties of observed 
massive early- type galaxies. Specifically, we simulate a 
Mr — 5\ogh = —22.5 galaxy with a half-light radius of 
12h~^ kpc and Sersic parameter of n = 4. Galaxies were 
assigned colors assuming a passively evolving simple stel- 
lar population (SSP) that was formed in a single burst 
at z = 3. For each galaxy, we convolve the simulated 
postage stamp with the local seeing, apply the flat field, 
bias, and bad column corrections in reverse, and add 
it to a raw SDSS image. Each image is then reduced 
using the standard SDSS PHOTO pipeline. Figure 1 
shows the result of this test. We find no significant trend 
with redshift of the measured flux compared to the total 
galaxy flux, indicating that our photometry is not biased 
strongly due to sky subtraction errors. The mean flux ra- 
tio found in our simulations, 80%, is quite close to that 
expected as the Petrosian flux systematically estimates 
the total flux of a galaxy with a n = 4 surface bright- 
ness profile to be ~ 82% of its total flux (Graham et al. 
2005). Throughout this work, we use the luminosity de- 
rived from the measured Petrosian fltix directly, and thus 
if comparisons are done to luminosity functions based on 
total flux measurements, care must be taken to account 
for this systematic effect. 
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Fig. 1. — Ratio of reconstructed Petrosian flux to the total galaxy 
light for 2,000 simulated galaxies with Mr — 5\ogh = —22.5, half- 
light radii of 12/i^^kpc, and colors of a passively fading SSP formed 
at z = 3. The dark asterisks mark the mean and la dispersion of 
the simulations while the gray points show each of the fake galaxy 
trials. We find no mean trend in the recovered flux with redshift 
and thus our galaxies are unaffected by overestimates of the local 
sky background which lead to underestimated galaxy fluxes for 
very large galaxies at low redshift. 

While Petrosian fluxes are unbiased in the absence of 
seeing, as a galaxy becomes unresolved, the Petrosian 
flux will report a systematically smaller fraction of the 
galaxy light (Blanton et al. 2001). Similarly, when work- 
ing near the detection limit of our imaging, one may 
worry that a given object only scatters above the de- 
tection threshold a fraction of the time; an average flux 
across many epochs can systematically overestimate the 
flux of such a source. At 2; > 0.7, the sizes of our sam- 
ple galaxies are approaching the size of the typical SDSS 
seeing disk and are quite faint relative to typical SDSS 
applications. To ensure that photometry of these high- 
redshift galaxies are unbiased, we simulate 10,000 galax- 
ies dX z > 0.7 with Mr — 51ogft, = —21.5 (corresponding 
to the faintest galaxies used in our luminosity function 
calculations in §3), half-light radii of 8/i~^ kpc, and col- 
ors characteristic of a passively fading SSP which formed 
at z = 3. Using an identical procedure to that described 
in §2.2, we add simulated images to raw SDSS frames 
and measure their photometry using PHOTO. We gen- 
erate 30 realizations of the simulations with the galaxy 
parameters and positions held constant but allowing the 
Poisson noise of the fake stamp to vary between realiza- 
tions. We tlic^n coadd the photometric measurements in 
each fake observation epoch to generate a mock coadded 
catalog of massive high-redshift galaxies using the same 
method described in §2.1 to generate the SDSS coadded 
catalog. Figure 2 shows the results of this test for the 
SDSS z-band which is the basis of our high redshift lumi- 
nosity measurements. The grey points show each galaxy 
simulated in this experiment while the stars show the 
mean in bins of input total flux. The mean ratio of Pet- 
rosian flux to input total flux is consistent with the ratio 
of 80% measured for low-redshift simulations above and 
thus we do not expect our use of Petrosian quantities 
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Fig. 2. — Simulation of coadded Petrosian flux measurements 
in high-redshift photometric data. Each grey point represents 
the coadded Petrosian flux from 30 realizations measured with 
the same method used to coadd the individual SDSS photomet- 
ric epochs to generate our deep photometric catalog. The mean in 
bins of total flux are shown as stars. Each photometric galaxy has 
properties of known high-redshift massive galaxies and thus the 
input flux, color, and size are all correlated - the faintest galaxies 
in this figure arc also the smallest. Wc find that galaxies above 
the 2-band flux limit (vertical dashed line) arc not strongly af- 
fected by the seeing disk; the (;,r,and i bands follow similar trends. 
The horizontal dashed line shows the mean flux ratio measured for 
low-redshift simulations. 

when measuring luminosities to bias our results to the 
flux limit of our survey (shown by the vertical dashed 
line). Below our selection limit, galaxies become unre- 
solved and the total recovered flux begins to decline. 

2.3. High-redshift Galaxy Sample 

The 54s exposure time of SDSS imaging is not suffi- 
cient to select galaxies at 2; ~ 0.9 based on their colors. 

The added depth of the SDSS Southern Survey, however, 
allows for the selection of massive galaxies to z ~ 1.0. Us- 
ing a similar method utilized to select LRGs at moderate 
redshifts from SDSS, we employ color cuts in griz to iso- 
late high-redshift LRGs for spectroscopy. In designing 
this selection, we capitalize on the fact that the strong 
4OOOA break of early-type galaxies moves through the 
?'-band at 0.6 < z < 1 resulting in progressively redder 
i — z colors while the r — i color shows less variation. 
Figure 3 illustrates the expected color evolution of mas- 
sive galaxies at z > 0.5. The gray scale shows the locus 
of galaxy colors from the deep SDSS imaging. The solid 
curves show the expcictcid evolutionary tracks for three 
different star formation histories; the reddest curve in 
r — i is a very early-type SED while the bluest track in 
r — i is roughly an early-type spiral (e.g. an Sa) from 
Bruzual & Chariot (2003). Galaxies with later spectral 
types never get comparably red in r — i; for comparison, 
the dot-dashed track shows the color evolution of an Sc 
type galaxy. The open circles are separated by Az = 0.1 
with the break in the color tracks occurring at z ^ 0.7. 
Above z ~ 0.7, the r — z color measures the distance 
from the turn in the color tracks and thus provides a 
good estimate of the photometric redshift of early-type 



galaxies. 

We construct two regions in this color-color space to 
select galaxies for deep spectroscopic observations. Sim- 
ilarly to Eisenstein et al. (2001), we define 

C_L = (r - z)model - {9 - f')model/4 - 0.177. (2) 

We require every galaxy candidate to satisfy 

^psf ^modcl > 0.2, (3) 

0.15 <ci^< 1.2 (4) 

< (r - i)modc\ < 1.7 (5) 

0.3 <{t- z)^odei < 1.5 (6) 

17 < ^model < 20.3 (7) 

1.5 < (r - z)n,ode\ < 2.5 (8) 

Here, the magnitude and color subscripts mark if the 
magnitude was based on SDSS PSF magnitudes or 
MODEL magnitudes (Stoughton et al. 2002). Equation 
(3) limits targets to objects in which at least 20% of the 
flux arises outside a central point source to select only 
extended objects in the SDSS photometry. At z = 0.9, 
1.2 arcseconds (the median seeing of our deep photome- 
try) corresponds to 6.7 kpc, smaller than the typical 
luminous red galaxy, and thus we do not expect galaxies 
of interest to be unresolved at 0.7 < z < 0.9. The defini- 
tion of c_L follows that of Eisenstein et al. (2001) and is 
designed to be parallel to the low-redshift galaxy locus in 
g — r versus r — i color-color space; Equation (4) removes 
z < 0.45 interlopers from the sample. 
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Fig. 3. — Selection of massive red galaxies at z > 0.5. The greyscale 
illustrates the observed galaxy locus for galaxies brighter than ^modei = 
20.3 from the SDSS Southern Survey. The three solid tracks show the 
expected colors of passively fading galaxies from Bruzual & Chariot 
(2003). The reddest track in r — i shows the expected colors of a very 
early-type galaxy and the bluest solid track shows those of an early- 
type disk galaxy (such as an Sa). The dot-dashed track shows the 
colors of an Sc type galaxy, for comparison. The tracks arc marked 
by open circles at Az — 0.1 intervals between redshifts of 0.5 and 1.0; 
the strong break in the colors occurs at z ^ 0.7. The boxed regions 
illustrate our photometric color selection. As detailed in section §2.3, 
galaxies at i — z > 0.6 are targeted at higher priority than galaxies 
with 0.3 < i — z < 0.6 as the redder galaxies are most likely to reside 
a.t z > 0.7. 
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Fig. 4. — Angular correlation functions for stars and galaxies 
selected with our liigli-rcdshift galaxy color criteria. The galaxy- 
galaxy correlation function (squares) shows strong clustering on 
all scales while both the star-star auto-correlation function (aster- 
isks) and star-galaxy cross-correlation function (diamonds) show 
very little clustering signal on several arcminute scales. If many 
galaxies were lost from our sample due to being unresolved by our 
star-galaxy separation, the star-galaxy cross correlation function 
would mirror that of the galaxy-galaxy auto-correlation function. 
Thus, the lack of signal at small separations in the stax-galaxy cross 
correlation function indicates we lose, at most, 2% of our galaxy 
targets due to our star-galaxy separation errors. 

Equations (4) - (8) limit our sample to red galaxies at 
0.5 < 2: < 1.0 and the flux limit imposed by Equation (7) 
isolates only the most luminous galaxies in this redshift 
range. We divide our selection into two groups based on 
their i — z color. Galaxies with i — z > 0.6 are given higher 
priority than galaxies with 0.3 < i — z < 0.6 as the redder 
subset of galaxies are more likely to lie at z > 0.7 as 
shown in Figure 3. Based on early observations and data 
simulations, we found that our redshift success would 
degrade at fluxes fainter than Zmodci = 20. In order to 
maximize the number of high-quality rcdshifts obtained, 
we targeted galaxies at -Zmodei < 20 at a higher priority 
than galaxies with 20 < z^^dei < 20.3. After target 
selection, fibers were allocated to 20% of the available 
galaxy candidates in the field. 

If there are unresolved galaxies that were untargetted 
with our algorithm, we can quantify this sample bias by 
comparing the galaxy angular correlation function to the 
star-galaxy cross correlation function from our targeting 
data. As the locations of distant galaxies are uncor- 
related with Galactic stars, the presence of unresolved 
galaxies in our star sample will result in an apparent 
signal in the star-galaxy cross-correlation function due 
to the correlated galaxy interlopers in the sample. We 
construct a sample of stars which meet identical selec- 
tion criteria used to select galaxies with the exception 
of Equation (3). After masking out 2' regions around 
bright (r < 12) stars, we count the number of galaxy- 
galaxy, star-galaxy, and star-star pairs as a function of 
separation compared to the expected number of pairs de- 
rived from a mock catalog of objects over the same area 
and subject to the same bright star mask. Our spectro- 
scopic observations directly probe the contamination by 



stars in our galaxy sample; we use this known contamina- 
tion rate to correct for the dilution of the galaxy-galaxy 
auto-correlation function arising from the addition of an 
uncorrelated stellar sample and create the average cor- 
relation function shown in Figure 4. As expected, the 
star-star auto-correlation function (asterisks) shows lit- 
tle power on several arcminute scales whereas the galaxy- 
galaxy auto-correlation (squares) function shows signif- 
icant clustering. The lack of strong signal in the star- 
galaxy cross-correlation function implies only a small 
fraction of galaxies can be lost to the star sample. Based 
on our measurements, we find that a maximum of 3% of 
the star sample can be contributed by interloper galax- 
ies at 99% confidence. As the average number density of 
stars in our fields is about 40% larger than galaxy tar- 
gets, we find that we lose, at most, 2% of our galaxy 
targets due to our star-galaxy separation. 

2.4. MMT Spectroscopy Observations and Data 

Processing 

We observed selected galaxies using Hectospec (Fab- 
ricant et al. 1998, 2005; Roll et al. 1998), a 300-fiber 
spectrograph on the 6.5m MMT telescope between Mar 
2004 and Oct 2005. Hectospec oflFers a 1 deg^ field of view 
and covers from 4000-9000 A with 6A resolution. Obser- 
vations were completed tising seven pointings with Hec- 
tospec. For each field, approximately half of the fibers 
were used to target high-redshift massive red galaxy can- 
didates and half were used to measure the faint quasar 
luminosity function (Jiang et al. 2006). Exposure times 
varied due to conditions, but each field was observed for 
an average of 3 hours. 
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Fig. 5. — Example of MMT spectra of high-redshift galaxies. 
Each spectrum has been smoothed by 2 resolution elements for dis- 
play; the spectra each have resolution of 6A. In each panel, vertical 
lines highlight prominent spectral features to guide the eye. The 
top panel shows a z = 0.92 galaxy with moderate signal-to-noise. 
The strong Call H-l-K absorptions lines and G-band at 4300A al- 
low for accurate redshift determination even at low signal-to-noise 
ratio. The middle panel shows a high signal-to-noise ratio z = 0.76 
spectrum and the bottom panels shows a 2 = 0.82 galaxy with 
strong Balmer absorption features characteristic of 1 Gyr popula- 
tions. The spectral range plotted was chosen to highlight the key 
features of our spectra; Hectospec observes considerably further 
into the blue but those data arc generally of quite low signal-to- 
noise for the high-redshift galaxies studied here. 
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All Hectospec data were reduced using the HSRED 
^ package which is based upon the SDSS spectroscopic 
pipeline. Data were flat-fielded using observations of an 
illuminated screen in the dome to remove pixel-to-pixel 
sensitivity variations as well as to correct for the strong 
fringing in the Hectospec CCDs in the red. When possi- 
ble, spectra of the twilight sky were taken to provide a 
secondary correction to account for any low-order resid- 
uals between fibers after the flat field derived from the 
dome flat corrections were applied. Wavelength solutions 
were obtained each night using observations of HeNeAr 
calibration lamps and the location of strong emission 
lines in the spectrum of the night sky were used to correct 
for any drift in the wavelength solution between the ob- 
servations of the calibration frames and the data frames. 

Observations of each field included approximately 30 
sky fibers which we used to construct the master sky 
spectrum from each exposure and subtract that from 
each object spectrum. Additionally, 3-5 photometrically 
selected F stars were targeted in each field. The ex- 
tracted spectra of these stars are compared to a grid 
of Kurucz (1993) model atmospheres to determine the 
spectral type of each star. Once we have determined 
the spectral type of each F star, we measure the aver- 
age ratio between the observed spectra and the model 
prediction to determine the global calibration to convert 
counts pixel" ^ to ergs s~^ cm~^ A~^. Figure 5 shows 
three fully-processed spectra from this survey. 

To determine the redshift of each object we compare 
the observed spectra with stellar, galaxy, and quasar 
template spectra and choose the template and redshift 
which minimizes the between model and data. As 
many of our spectra have low signal-to-noise ratios, ev- 
ery spectrum is examined by eye to ensure that the fitted 
redshift was correct. In cases in which the automated 
routine failed to converge to the correct redshift, a hand- 
measured redshift is used in its place. Our spectroscopy 
resulted in redshifts for 470 galaxies at 0.6 < z < 1.0 
over 7 deg^ and 302 galaxies at 0.7 < z < 0.9 which will 
be used in our analysis, here. Figure 6 shows the color 
distribution of the confirmed galaxies at 0.7 < z < 0.9 
which are used for our luminosity function calculations 
at high redshift. Of the 890 galaxy candidates that were 
targeted for spectroscopy, 12% of the spectra did not re- 
sult in a redshift measurement. 

3. LUMINOSITY FUNCTION CONSTRUCTION 

3.1. Calculation of Rest-frame Luminosities 

In order to compare the populations of massive red 
galaxies as a function of redshift, we first need to trans- 
form the observed photometry to the rest-frame of each 
galaxy to remove the effects of redshift on the observed 
properties. A number of approaches have been developed 
to perform fc-corrections to the rest-frame system; each 
approach has its advantages and drawbacks. In order 
to minimize errors introduced due to errors in the stel- 
lar synthesis models used to calculate our fc-corrections, 
we consider the rest-frame properties of our galaxies 
through a modified SDSS filter set. This system, de- 
noted 0.3^0. 3^0. 3^0. 3^0.3^^ consists of the SDSS ugriz fil- 
ters which have been blueshiftcd by a redshift of 0.3 sim- 
ilar to the approach used in Blanton et al. (2003), Cool 

^ http:/ /mizar.as.arizona.cdu/hsrcd/indcx.html 




0.8 1.0 1.2 1.4 1.6 
Observed r-i 

Fig. 6. — Colors of confirmed galaxies at 0.7 < z < 0.9 from 
our MMT spectroscopy. The early-type galaxy color tracks and 
color selection criteria are as shown in Figure 3. The colored 
points show the location of each of our sample galaxies in this 
color space; the shape (color) of each point denotes its redshift. 
Stars (magenta) show 0.70 < z < 0.75 galaxies, diamonds (blue) 
mark 0.75 < z < 0.80 objects, and the squares (green) and filled 
circles (red) illustrate 0.80 < z < 0.85 and 0.85 < z < 0.90 galax- 
ies respectively. The grey dot-dashed lines show the sub-regions of 
color-color space used to measure the fraction of spectroscopically 
observed galaxies which were excluded when evolved to our lower- 
redshift bins. We use this correction factor when bootstrapping to 
our full photometric sample as described in §3.2. 

et al. (2006), and Wake et al. (2006). In this system, 
a galaxy at a z = 0.3 will have a fc-correction that is 
independent of its spectral energy distribution and will 
equal — 2.51ogiQ(l-|-0.3). We choose a shift of 0.3 to draw 
upon the fact that at z ~ 0.8 (near the median redshift 
of our high- redshift galaxy sample), the observed z-band 
probes a similar portion of the spectrum as probed by 
the r-band observing a z = 0.3 galaxy. In the following 
sections, we will measure the Mo. a,, luminosity function 
of massive galaxies; for comparison, B —^-^ r « —0.01 for 
an old stellar population. Based on luminosity function 
fits from Brown et al. (2007), Ml^^ - 51og/i = -20.3 and 
thus our sample focuses on galaxies with L > 3L*. For 
reference, a 3L* SSP at z = 0.3 which formed its stars 
at z = 3 has an approximate stellar mass of 3 x 10^^ Mq. 

To construct the fc-corrections for galaxies in each 
of our samples, we create a grid of evolving and non- 
evolving SSP at solar metallicity with formation redshifts 
ranging from 1 to 10 from Bruzual & Chariot (2003) 
based on a Salpeter (1955) initial mass function (IMF). 
We find that this set of models adequately span the range 
of observed colors for all of our galaxies. Each galaxy 
is assigned a template based on a maximum likelihood 
comparison of the predicted colors and observed SDSS 
photometry. 

While the ^-corrections based on non-evolving models 
assume that the underlying stellar population remains 
unchanged from the observed epoch, our k-\-e corrections 
include the passive evolution, normalized to z = 0.3, 
of the stellar populations in the galaxies between the 
observed epoch and the rest-frame redshift. For each 
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FlK. 7. Rcdshift success versus the z-band flux in a 1.5 arc- 
second aperture for two of our targeted fields. The triangles show 
a high quality mask observed under photometric conditions and 
excellent (ft; 0.5") seeing. The asterisks show a poor-quality mask 
affected by clouds and poor seeing leading to degraded success 
at the faintest fluxes. We correct for this incompleteness in each 
of our Hectospcc fields before computing the luminosity function 
using low-order fits as show by the dot-dashed lines. The grey his- 
togram illustrates the distribution of fiber magnitudes for all of our 
spectroscopic targets. The sharp decline in objects at ^flber = 20.8 
corresponds to our sparser sampling of objects with Zmodel > 20. 
The squares mark the redshift completeness of our full spectro- 
scopic sample. 

galaxy, wc use the best fitting SSP to predict the SED 
the galaxy would have a,t z — 0.3; a galaxy fit by a SSP 
with age t will age into a SSP with age t -I- At{zo) where 
At{zo) is the lookback time difference between z = 0.3 
and zo, the observed redshift of the galaxy. We include 
both types of models in order to compare the affects of 
passive evolution on the inferred evolution of the lumi- 
nosity function of massive galaxies since z ~ 0.9. 

3.2. Luminosity Functions 

Luminosity functions are calculated using the standard 
method (Schmidt 1968). For each galaxy, we cal- 
culate the redshifts at which the galaxy would have been 
selected and observed in our survey. In this calculation, 
we utilize the best-fit template chosen when calculat- 
ing fc-corrcctions, as described above, to estimate each 
galaxy's colors as a function of redshift. Based on these 
predicted colors, we assign a probability (0 or 1) that a 
given galaxy would have been selected at each redshift. 
The maximum available volume is then the integral over 
the redshift range weighted by the selection probability 
at each redshift. 

Each sample is corrected independently for the spectro- 
scopic completeness of the observations. The low-redshift 
SDSS MAIN and intermediate-redshift SDSS LRG galax- 
ies were corrected to account for the spatially-dependent 
incompleteness of SDSS spectroscopy. As we have sev- 
eral priority classes in our high-redshift target selection, 
we must correct our sample with more detail than merely 
the fraction of the galaxies that received fibers. Instead, 
we break our sample into four regions in color- magnitude 
space and calculate the completeness in each region in- 



dependently. As described in §2.2, galaxies were given 
priority based both on their i — z color and Zmodoi flux. 
This results in four color-magnitude regions in which we 
then calculate the photometric completeness by counting 
the number of photometrically selected galaxies which 
were given a fiber compared to the number of galaxies in 
the parent catalog in that color and magnitude bin. Our 
completeness correction was calculated independently for 
each of our seven Hectospec fields. In each field, we com- 
pare the number of spectroscopically observed objects to 
the total number of photometric objects within a 2 deg^ 
square box around the field center when calculating our 
incompleteness. In doing this, we bootstrap our spectro- 
scopic sample to 9000 photometrically selected galaxies 
over twice the area observed with Hectospec, thus mini- 
mizing the effects of cosmic variance on our sample. The 
inclusion of this photometric sample doesn't change the 
normalization of the high-redshift luminosity function we 
measure, but results in smaller errors due to field-to-field 
variations in the galaxy number counts. 

Signal-to-noisc ratio variations in our high-redshift 
galaxy spectroscopy result in approximately 12% of our 
observed objects with no measurable rcdshift. In order 
to correct for this effect, we measure the fraction of ob- 
served galaxies with viable redshifts as a function of the 
z-band flux within an 1.5 arcsecond aperture centered 
on our fiber location to estimate the flux available to the 
spectroscopic fiber. We then fit this relationship with a 
low-order polynomial for each Hectospec field and apply 
the derived correction before calculating the the final lu- 
minosity function. Figure 7 shows an example of this 
technique on two different fields spanning the full range 
of data quality. The triangle symbols show the com- 
pleteness for a field with high signal-to-noise observed 
under photometric conditions and superb seeing O'.'S) 
while the asterisks show a field observed under less pho- 
tometric conditions. The range in data quality leads to 
significant c;ompletcness variations between each of our 
spectroscopically observed fields; neglecting this would 
bias our final inferred luminosity function. The square 
symbols in the figure show the composite completeness 
for the full galaxy sample as a function of fiber magni- 
tude. 

We make a further correction to ensure that the galax- 
ies utilized in the construction of the luminosity function 
in each redshift bin probe a homogeneous population of 
objects. Using the best-fit stellar population template 
derived when calculating the k + e-corrections, wc es- 
timate the colors of each galaxy as a function of red- 
shift from z = 0.1 to z = 0.9. We then require that ev- 
ery galaxy included in our calculation of the luminosity 
function would have been selected in each of our redshift 
samples thus ensuring that the population of galaxies wc 
consider at 0.1 < z < 0.2 arc consistent with galaxies at 
0.7 < z < 0.9 after the passive evolution of their stel- 
lar populations has been included. When bootstrapping 
to the entire photometric sample of galaxies at high red- 
shift, we grid the r—i versus i — z color-color plane into 12 
subsections as shown in Figure 6 and calculate the frac- 
tion of galaxies in each subrcgion that would be excluded 
based on this criterion. The size of these sub regions was 
chosen to sample both the i — 2: < 0.6 and the i — z > 0.6 
subsamplcs with similar detail. The final results are not 
strongly dependent on the exact subregions chosen for 
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this correction. 

When target selection is based on noisy photometry, 
the effects of photometric scattering of objects into or 
out of the nominal color- and flux-limits can be quite 
significant (Wake et al. 2006). As our high-redshift sam- 
ple of galaxies is selected from SDSS stacked photometry, 
we perform an empirical test of this photometric scatter- 
ing on our sample. Using the full sample of SDSS main 
galaxies observed at 0.1 < z < 0.2 we create a mock 
sample of 0.7 < z < 0.9 galaxies based on the best-fit 
k + e-corrections described in §3.1. We then subject this 
mock galaxy sample to representative photometric errors 
present in our coadded photometric catalog and deter- 
mine the fraction of mock galaxies that would have been 
selected in the presence of photometric errors. For galax- 
ies brighter than z = 20, we find that ~ 2% of selected 
galaxies have colors that would fall outside our color-cuts 
but scatter into the sample when photometric errors are 
included. At fainter magnitudes, 20 < z < 20.3, approx- 
imately 10% of the galaxies included in the mock high- 
redshift galaxy sample have scattered above the survey 
flux- limit due to photometric errors. When calculating 
our high-redshift luminosity functions, we include these 
contamination rates as a statistical weight assigned to 
each galaxy based on its observed z-band flux. 

In order to estimate the error on our high-redshift lu- 
minosity function measurements, we remove each of our 
spectroscopic fields (and ancillary photometric data), in 
turn, from our calculation of the ensemble luminosity 
function and repeat our calculations; we use the mea- 
sured variation in the luminosity functions created with 
this test as an estimate of the large scale structure error 
on our luminosity function measurements. Similarly, for 
our SDSS samples, we divide the SDSS survey area into 
20 subregions and perform the same experiment. These 
jack-knife errors are ~ 25% larger than those based on 
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Fig. 8. — Luminosity function of massive galaxies with only a 
fc-correction applied to account for the redshifting of galaxy light. 
The symbols (color) mark the four redshift bins used ; diamonds 
(black) 0.1 < z < 0.2, asterisks (magenta) 0.2 < 2 < 0.3, squares 
(green) 0.3 < z < 0.4, and circles (red) 0.7 < 2 < 0.9. The 
luminosity functions show the characteristic brightening toward 
higher redshifts due to the passive aging of stars. We must correct 
for the passive evolution of stellar populations in order to measure 
the evolution in the underlying galaxy population. 



Poisson errors alone in the lowest-luminosity bins and are 
comparable to those estimated from counting statistics 
at the bright end. While subsampling can result in an 
underestimate of the error if a single large scale feature is 
present in multiple subfields, the large area surveyed by 
SDSS at low redshift and the several degree separation 
between our spectroscopic fields at high-redshift mini- 
mize this effect and thus jack-knife errors are a robust 
estimate of the cosmic variance errors for our samples. 
Throughout this paper, we utilize the larger of the two 
errors when doing calculations with our measured lumi- 
nosity functions. 

Figure 8 and Table 1 show the non-evolving lumi- 
nosity function measured from our samples. The sym- 
bol (color) denotes the redshift bin : diamonds (black) 
0.1 < z < 0.2, asterisks (magenta) 0.2 < z < 0.3, squares 
(green) 0.3 < z < 0.4, and circles (red) 0.7 < z < 0.9. 
The figure shows a clear separation between each lu- 
minosity function with higher-redshift galaxies having 
higher luminosities (or larger number density). This 
characteristic behavior is expected due to the passive fad- 
ing of the stellar populations in these massive red galax- 
ies. We must remove this effect in order to understand 
any true changes in the underlying population of massive 
galaxies since z ~ 0.9. The turnover at low-luminosities 
is an artifact of the color-selection of these galaxies. As 
shown in Eisenstein et al. (2001), the LRG sample selec- 
tion results in a diagonal cut across the red-sequence at 
low luminosities which is being refiected here as the turn 
over at low-luminosities in our luminosity function. This 
should not be interpreted as a characteristic luminosity 
of the sample. The luminosity functions of galaxies in 
our survey are shown in Figure 9 and recorded in Ta- 
ble 2 after the effects of evolution are included. After 
the effects of passive evolution are accounted for, the lu- 
minosity functions show little variation between redshift 
bins. The integrated luminosity densities for both the 
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Fig. 9. — Luminosity function of massive galaxies after both the 
redshifting of their spectra and the passive evolution of their stel- 
lar populations have been accounted for when calculating galaxy 
luminosities. The symbols are as described in Figure 8. We find 
very little evolution in the number counts of massive galaxies to 
2 ~ 0.9, indicating that the most massive galaxies have grown little 
over the latter half of cosmic history. 
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TABLE 1 

Lriiixors Red C-Ialax^" LriiL\A)siTA" l-'rxA iioxs \\"iTn A'o E\A)LrTioxAin" Couhfatiox 



logio Galaxy Number Density'' logio Galaxy Number Density'' 

Mo.3^-5logh 0.1<z<0.2 0.2<2<0.3 0.3 < z < 0.4 Mo.a^-Slog^i 0.7 < z < 0.9 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


-21.55 


-3.66 ±0.01 


-3.64 ± 0.03 


-3.78 ± 0.05 


-21.59 


-3.26 ± 0.08 


-21.65 


-3.75 ±0.01 


-3.56 ± 0.07 


-3.66 ±0.04 


-21.77 


-3.30 ±0.07 


-21.75 


-3.89 ±0.01 


-3.85 ± 0.01 


-3.64 ±0.04 


-21.95 


-3.44 ±0.06 


-21.85 


-4.04 ±0.02 


-3.92 ± 0.04 


-3.73 ±0.03 


-22.13 


-3.53 ±0.06 


-21.95 


-4.21 ±0.02 


-4.06 ± 0.03 


-3.87 ±0.02 


-22.31 


-3.89 ±0.09 


-22.05 


-4.38 ± 0.03 


-4.25 ± 0.01 


-4.02 ±0.02 


-22.49 


-4.28 ± 0.13 


-22.15 


-4.53 ±0.04 


-4.41 ± 0.02 


-4.18 ±0.01 


-22.67 


-4.64 ± 0.24 


-22.25 


-4.79 ±0.04 


-4.57 ±0.02 


-4.34 ±0.02 


-22.85 


-5.40 ± 0.43 


-22.35 


-4.94 ±0.05 


-4.77 ±0.02 


-4.63 ±0.01 






-22.45 


-5.15 ±0.06 


-4.97 ± 0.03 


-4.83 ±0.01 






-22.55 


-5.36 ±0.08 


-5.16 ± 0.04 


-5.01 ±0.02 






-22.65 


-5.65 ±0.12 


-5.38 ±0.05 


-5.19 ±0.02 






-22.75 


-5.74 ±0.13 


-5.78 ±0.07 


-5.41 ± 0.03 






-22.85 


-6.18 ±0.22 


-5.88 ± 0.08 


-5.68 ± 0.03 






-22.95 


-6.13 ±0.22 


-5.99 ±0.09 


-5.88 ± 0.04 







"■All number densities are expressed in units of Mpc ^ Mag 



TABLE 2 

Luminous Red Galaxy Luminosity Functions After Passive Evolution Correction 



logio Galaxy Number Density'' logio Galaxy Number Density" 



Mo.3r-5 1og/i 0.1 < 2 < 0.2 0.2<2<0.3 0.3 < 2 < 0.4 Mo.3r--5logh 0.7 < 2 < 0.9 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


-21.55 


-3.52 ±0.01 


-3.64 ±0.07 


-3.71 ±0.04 


-21.59 


-3.48 ± 0.06 


-21.65 


-3.56 ±0.01 


-3.60 ± 0.02 


-3.62 ±0.04 


-21.77 


-3.68 ± 0.07 


-21.75 


-3.66 ±0.01 


-3.60 ±0.03 


-3.66 ±0.04 


-21.95 


-3.96 ± 0.11 


-21.85 


-3.74 ±0.01 


-3.72 ± 0.02 


-3.81 ±0.02 


-22.13 


-4.25 ± 0.12 


-21.95 


-3.90 ±0.01 


-3.90 ±0.01 


-3.94 ±0.02 


-22.31 


-4.55 ± 0.16 


-22.05 


-4.05 ±0.02 


-4.05 ± 0.01 


-4.12 ±0.01 


-22.49 


-4.79 ± 0.22 


-22.15 


-4.21 ±0.02 


-4.21 ± 0.01 


-4.24 ±0.02 


-22.67 


-5.41 ± 0.43 


-22.25 


-4.40 ± 0.03 


-4.37 ±0.02 


-4.45 ±0.01 


-22.85 


-5.40 ± 0.43 


-22.35 


-4.57 ± 0.03 


-4.57 ±0.02 


-4.61 ±0.02 






-22.45 


-4.75 ±0.05 


-4.74 ± 0.03 


-4.81 ± 0.02 






-22.55 


-5.06 ±0.06 


-4.92 ± 0.03 


-4.97 ±0.03 






-22.65 


-5.23 ±0.07 


-5.12 ±0.04 


-5.17 ±0.03 






-22.75 


-5.33 ±0.08 


-5.50 ± 0.06 


-5.42 ± 0.04 






-22.85 


-5.73 ±0.13 


-5.70 ± 0.08 


-5.82 ± 0.07 






-22.95 


-5.89 ±0.15 


-5.85 ± 0.09 


-5.96 ± 0.08 







"All number densities arc expressed in units of h'^ Mpc Mag 



evolutionary-corrected and fc-corrected luminosity func- 
tions are listed in Table 3. Analysis of these luminosity 
functions is the focus of §4. 

4. LUMINOSITY EUNCTION ANALYSIS 

4.1. Evolution in the Massive Galaxy Population Since 

z ~ 0.9 

The agreement between the luminosity function mea- 
surements at 0.1 < z < 0.9 as illustrated in Figure 9 in- 
dicates that the massive galaxy population has evolved 
little since z ~ 0.9. In order to quantify this evolution, 
we have adopted a similar parameterization to that dis- 
cussed by Brown et al. (2007). Instead of measuring 
the evolution in the total luminosity density contained 
in massive galaxies, we instead measure the magnitude 
at which the integrated number density reaches a certain 
value. As massive galaxies populate the exponential tail 
of the luminosity distribution, small photometric errors 
can result in significant errors in the total luminosity 



density derived. For example, a shift of 3% in the lu- 
minosity threshold corresponds to a 10% change in the 
inferred number density of the population. Thus, if the 
integrated number or luminosity density at a given mag- 
nitude is used to measure the evolution of a population, 
results are quite sensitive to the magnitude threshold uti- 
lized. Here, we use the inverse; we measure the magni- 
tude at which the integrated number density reaches a 
threshold of 10-4-5 and 10-^ ° Mpc'^. These mag- 
nitudes are denoted by Mo.3r{10~'^-^) and Mo.3r{10~^-°) 
throughout this discussion. 

In order to measure Afo.3^(10^''"'') and Mo,3,.(10^'''°), 
we fit each of our luminosity functions with a quadratic 
polynomial in the logarithm. We then integrate the 
best fitting polynomial and determine the magnitude at 
which the integrated number density reaches 10-4-5 and 
10~^ '^h^ Mpc^"^. Error bars were calculated by repeat- 
ing this calculation while removing one of our subfields in 
turn in the same manner we calculated jack-knife errors 
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TABLE 3 

IXTECHATI' I) Lr.MIXOSIT^" DEXSIT^" 



j(Mo.3^ < -21.50)^ jjMo.s^ < -22.25)°- 



Rcdshift fcb 
(1) (2) 



0.15 
0.25 
0.35 
0.80 



2.63 ± 0.04 
3.54 ± 0.06 
5.18 ±0.11 
34.32 ± 10.00 



k + e" 
(3) 



5.45 ± 0.06 
5.22 ±0.09 
4.58 ± 0.09 
5.47 ±1.00 



(4) 



0.24 ±0.02 
0.33 ± 0.01 
0.50 ±0.01 
2.46 ± 0.62 



k + e" 
(5) 



0.54 ±0.03 
0.54 ±0.02 
0.48 ± 0.01 
0.71 ±0.29 



Luminosity densities in units of lO^h^ Lq Mpc ^ 

^ Integrated luminosity densities based on luminosity functions de- 
rived without correcting for the passive fading of stellar populations 

" Integrated luminosity densities based on luminosity functions calcu- 
lated after correcting for stellar evolution based on Bruzual & Chariot 
(2003) stellar population synthesis models and a Salpeter IMF as dis- 
cussed in §3.2. 



on our luminosity function measurements. The exact 
form we use to fit the luminosity function has little effect 
on our final results. Figure 10 shows the evolution in 
Mo-3r{lO^'^'^) and A'/o. 3^(10^^''') before the passive evo- 
lution of stellar populations is removed from our galaxies 
and columns (2) and (4) of Table 4 reports these mea- 
surements. Columns (2) and (6) of Table 5 and Figures 
11 and 12 show the same critical magnitudes recalculated 
after the affects of passive evolution have been removed 
from our galaxy luminosity measurements. In both fig- 
ures, the differences between the number density mea- 
sured in each redshift bin are significant within our er- 
rors. The large area probed by SDSS makes cosmic vari- 
ance between the redshift bins smaller than the observed 
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Fig. 10.— The evolution oi Ma.3^{10-'^-^) and Mo. 3^(10"^) , the 
magnitudes at which the integrated luminosity density reaches val- 
ues of 10~''-^/i'^ Mpc~^ (asterisks) and lO^^ '^/i'^ Mpc^-^ (squares) 
respectively. Here, we show the evolution of this parameter if the 
passive fading of stellar populations is not removed when calculat- 
ing galaxy luminosities. Both measurements show the character- 
istic brightening toward higher redshifts. Without removing the 
luminosity evolution induced by the passive evolution of stars in 
these massive galaxies, the observed trends may be due to both the 
passive fading of galaxies over time or the build up in the number 
density of these galaxies over cosmic history. 
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Fig. 11. — The evolution of Mo.3^(10~''-^), the magnitude at 
which the integrated luminosity function reaches a number den- 
sity of 10~* ®/i Mpc~^. This parameter is used to quantify the 
evolution of the LRG population as these galaxies populate the ex- 
ponential tail of the luminosity function and small changes to the 
magnitude threshold chosen may lead to significant errors when cal- 
culating the total number or luminosity density in these objects. 
The asterisks show measurements using the Bruzual & Chariot 
(2003) stellar templates, the squares show the derived evolution 
based on Maraston (2005) models (see §4.2), and the diamonds 
show measurements based on the flux within fixed 20h~^ kpc aper- 
tures and Bruzual & Chariot (2003) k + e corrections as described 
in §4.5. For clarity, the Bruzual & Chariot (2003) and Maraston 
(2005) points have been shifted by -0.02 and +0.02 in rcdshift, 
respectively. None of these samples shows a strong evolution in 
the massive galaxy population since z = 0.9. The dot-dashed line 
shows the best fit linear relationship based upon the Bruzual & 
Chariot (2003)-derived luminosity functions and the shaded area 
shows the l-a confidence of the fit. The best fitting slope predicts 
an evolution of 0.03 ± 0.08 mag between 2 = and 2 = 1 and is 
consistent with no-evolution (shown by the dotted line). 

differences at 0.1 < z < 0.4, so large scale structure is 
unlikely the cause. Wc fit the measured critical magni- 
tudes with a linear evolution with rcdshift. The best fit 
relation is shown as dot-dashed lines in Figures 11 and 
12; the shaded region shows the one sigma confidence of 
the fit. Fits to both critical magnitude thresholds find 
similar evolution; the critical magnitudes have evolved 
by 0.03 ± 0.08 mag between z = and z ~ 1. When 
fitting this value, we add systematic floor of 0.02 mag in 
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quadrature to each magntiude threshold. As shown by 
the dotted lines in the figures, the best fit to our data 
does not rule out pure passive evolution in the massive 
galaxy population, 
a 

4.2. Importance of k-corrections on the Result 

Central to any study of the rest-frame photometric 
properties of extragalactic sources are the fc-corrections 
used to convert the observed quantities to the rest-frame 
properties of the galaxy. There are a number of inherent 
problems with this method, in particular when applied 
to the massive galaxies of interest here. As demonstrated 
in Eisenstein et al. (2003) and Cool et al. (2006), pop- 
ular stellar synthesis models such as Bruzual & Chariot 
(2003) and PEGASE.2 (Fioc & Rocca-Volmerange 1999) 
do not match the spectral properties of LRGs, especially 
a-element features; LRGs are a-enhanced compared to 
solar while the synthesis models do not include non-solar 
a-abundances. Furthermore, a number of studies (e.g. 
Eisenstein et al. 2001; Wake et al. 2006) demonstrate that 
the current generation of stellar synthesis models poorly 
reconstruct the observed broad-band colors of galaxies 
on the red-sequence over a variety of redshifts. 

To explore the importance of the fc-correction mod- 
els on our inferred results, we employ a second set of 
fc-corrections based on the Maraston (2005) models pro- 
vided by C. Maraston (private communication). These 
models were created to more accurately track the colors 
of massive red galaxies than simple stellar populations. 
The spectrum is modeled as a composite of a metal-rich 
(2Zq) population and a metal poor (0.005Zq) popula- 
tion; the metal-poor population holds 10% of the mass 
in the galaxy. 

Figure 13 shows the expected colors of a passively fad- 
ing galaxy from the the Maraston (2005) and Bruzual & 
Chariot (2003) models utilized in our analysis. As shown 
in the figure, at z > 0.6, the Maraston models predict 
significantly bluer g — r colors, and more closely follows 
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Fig. 12. — Same as Figure 11 except showing the evolution 
of Mo.3^(10~^"). The best fit to the k + e-corrected luminosity 
functions based on Bruzual & Chariot (2003) models is shown, 
again. The fit here is independently calculated from the one in 
Figure 11, but shows the same slope. 
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Fig. 13. — Predicted passively evolving color tracks from Bruzual 
& Chariot (2003) (solid line) and a composite stellar population 
based on Maraston (2005) models (dot-dashed) as described in §4.2. 
The data show the colors of galaxies in our intermediate and high- 
redshift samples. The Maraston (2005) models predict significantly 
bluer g — r colors at high redshifts which follow the observed locus 
of galaxy colors more closely than Bruzual & Chariot (2003) SSP 
predictions. 

the observed color locus of galaxies in our sample. While 
the g — r and g — i colors of galaxies are better matched 
with the Maraston (2005) models, the r — i colors pre- 
dicted from both templates are systematically bluer than 
observed galaxies. 

In order to understand any systematics introduced 
based on the stellar synthesis models used, we re- 
performed our analysis using the Maraston (2005) mod- 
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Fig. 14. — Evolution of the massive galaxy luminosity function 
using Maraston (2005) models when correcting for the redshift- 
ing of the galaxy spectra and the passive evolution of their stellar 
populations. The data points are as in Figure 8. The dot-dashed 
line shows the 0.1 < 2 < 0.2 luminosity function calculated using 
Bruzual & Chariot (2003) templates for comparison. We find no 
strong difference in the inferred evolution of massive galaxies when 
different stellar synthesis models are used. 
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TABLE 4 

Evolution of the Massive Red Galaxy Population Without Correcting for Passive 

Evolution 





Afo. 3^(10- 


-5-0) - 5 log/i=>- 


Afo. 3^(10- 


- 5 logfet- 


Redshift 


Petrosian 


20h~^ kpc Aperture 


Petrosian 


20h~^ kpc Aperture 




Luminosity 


Luminosity'^ 


Luminosity 


Luminosity'^ 


(1) 


(2) 


(3) 


(4) 


(5) 


0.15 


-22.08 ± 0.008 


-21.84 ±0.005 


-21.79 ±0.005 


-21.63 ±0.005 


0.25 


-22.16 ±0.010 


-21.91 ±0.010 


-21.87 ±0.010 


-21.69 ±0.010 


0.35 


-22.25 ±0.010 


-22.04 ± 0.010 


-21.98 ± 0.010 


-21.84 ±0.010 


0.80 


-22.60 ± 0.090 


-22.37 ±0.020 


-22.33 ± 0.052 


-22.20 ±0.033 



"Tlic magnitude at wliicli tlic integrated number density of LRGs reaelies 10 ''■'^h ^ 
''The magnitude at which the integrated number density of LRGs reaches 10"^'^h~"^ 



■^See §4.5 for a full description of the aperture luminosity functions 



Mpc 
Mpc 



TABLE 5 

Evolution of the Massive Red Galaxy Population After Correcting for Stellar Evolution 



Mo.3^(10-S-°)-5 1ogfe=' 



Redshift Petrosian 20h ^ kpc Aperture Maraston 
Luminosity Luminosity"^ Luminosity'^ 

(1) (2) (3) (4) 



Mo.3^(10-^-S) -5 logftb 



Petrosian 
Luminosity 

(5) 



20?i ' kpc Aperture 
Luminosity'^ 
(6) 



Maraston 
Luminosity'' 
(7) 



0.15 
0.25 
0.35 
0.80 



-22.27 ±0.008 
-22.30 ±0.010 
-22.25 ±0.010 
-22.28 ± 0.056 



-22.23 ± 0.007 
-22.26 ±0.010 
-22.24 ±0.010 
-22.25 ± 0.080 



-22.00 ±0.005 
-22.04 ±0.010 
-21.99 ±0.010 
-22.04 ± 0.060 



-21.99 ±0.005 
-22.02 ± 0.010 
-21.96 ±0.010 
-22.01 ± 0.053 



-21.81 ±0.005 
-21.85 ±0.010 
-21.80 ±0.010 
-21.81 ±0.037 



-21.94 ±0.005 
-21.97 ±0.010 
-21.94 ±0.010 
-21.96 ± 0.075 



''The magnitude at which the integrated number density of LRGs reaches 10~^ °fe~^ Mpc~^. 
''The magnitude at which the integrated number density of LRGs reaches 10~'^-^h~^ Mpc~^. 
"^See §4.5 for a description of the 20/i~'kpc aperture luminosity function. 

''Maraston luminosities are Petrosian flux measurements which have been k ± e-corrected using Maraston (2005) models and are 
described in §4.2 



els as the basis for our k— and k + e-corrections. Fig- 
ure 14 shows the result of this analysis compared to the 

low-rcdshift luminosity function derived using Bruzual 
& Chariot spectral templates. The number density of 
massive galaxies shows little evolution after the passive 
evolution of the stellar evolutions are taken into account 
regardless of the models used to perform the k + e- 
corrections as shown in Figures 11 and 12 and columns 
(4) and (7) in Table 5. There is, however, a net offset 
in the measured luminosity of galaxies between the two 
methods, so care must be taken that fc-correction differ- 
ences are accounted for when comparing galaxy samples 
from differing analysis techniques. To quantify any dif- 
ference in the implied evolution based on these two sets 
of stellar templates, we plot both the Bruzual & Chariot 
and Maraston derived Mo.3r{10~^-^) and Mo.3r(10"^ °) 
in Figures 11, 12. In both data sets, these quantities 
have only evolved by less than 0.05 mag since z ^ 0.9, 
implying that massive galaxies do little more than fade 
over the latter half of cosmic history. 

4.3. Merger Fraction from z 0.9 

Following the method dc^scribed in Wake et al. (2006), 
we construct a toy model for the merger history of LRGs 
to constrain the merger rate of massive red galaxies since 
2; ~ 0.9. Using our 0.1 < z < 0.2 luminosity function, we 
create a mock sample of galaxies and then allow a fixed 
fraction of them to have undergone a 1:1 merger since z = 
0.9. We then compare the luminosity function prediction 



for this mock sample to the observed luminosity function 
to determine the probability that both were drawn from 
the same population. Examples of predicted luminosity 
functions assuming different merger fractions are shown 
with the high-redshift data in Figure 15. 

Our high-redshift luminosity function is best fit by no 
merging over the latter half of cosmic history. Merger 
rates greater than 25% are ruled out with 50% confidence 
and merger rates larger than 40% are excluded at the 
99% level based on our measured high-redshift luminosity 
function. This result agrees with previous studies based 
on lower-redshift data and photometric redshift surveys 
(Brown et al. 2007; Masjedi et al. 2006, 2007; Wake et 
al. 2006). If less massive mergers are considered, more 
substantial merger rates are permitted. Performing the 
same experiment but instead considering 1:3 mergers, no 
merging is still favored, but rates as high as 40% are 
allowed at 50% confidence and only merger rates larger 
than 60% are ruled out at 99% confidence. These rate 
limits imply the total stellar mass in massive red galaxies 
from z ~ 0.9 must not have grown by more than 50% 
(at 99% confidence) in order to reproduce the observed 
luminosity functions. 

The fact that the most massive red galaxies appear to 
have evolved very little beyond the passive aging of their 
stellar populations since z ~ 0.9 is quite interesting. The 
most massive galaxies reside in the most massive dark 
matter halos - these halos have not remained static since 
z ~ 1. In a standard ACDM universe, the most massive 
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Fig. 15. — Models of the high-redshift luminosity function (points 
and errorbars) . Each of the soUd Unes shows a simulation in which 
our 0.1 < 2 < 0.2 luminosity function is evolved backward assum- 
ing a fixed fraction of the LRGs has doubled its luminosity through 
1:1 mergers between z ~ 0.9 and z ~ 0.1. Full details can be found 
in §4.3. Our data are consistent with no growth in the massive red 
galaxy population since z ~ 0.9 ; merger fractions larger than 25% 
are ruled out at the 50% confidence level and merger frac tions 
larger than 40% are ruled out at the 99% level. 

halos [M > 3 X VS^'^Mq) have grown by a factor of two 
or three since redshift of unity (Seo et al. 2007; Conroy et 
al. 2007b); one would naively estimate that the galaxies 
that reside in these halos would have grown, as well. 

LRGs at z = 0.3 are known to reside in dense environ- 
ments with mean clustering similar to rich groups and 
poor clusters (Zehavi et al. 2005). The formation and 
assembly of groups and clusters at z < 1 would naturally 
result in a discrepancy between the stellar mass growth 
of the massive central galaxy and the dark matter halo 
mass in which it resides. As satellite galaxies are accreted 
into the group or cluster halo, these satellites contribute 
stellar mass to the total stellar mass of the halo but not 
to the stellar mass of the central galaxy. The fact that 
galaxies with masses M > IO^^Mq are observed to re- 
side in a broad range of halo masses (Mcintosh et al. 
2007) may be a natural outcome of group and cluster 
formation. 

If the lack of evolution in the number density of LRGs 
is due to the growth of clusters rather than the growth 
of the central LRG, one would expect to observe multi- 
ple LRGs within a single cluster halo. To address this 
hypothesis. Ho et al. (2007) performed a thorough ac- 
counting of the number of LRGs which reside in a single 
halo in the SDSS dataset and Conroy et al. (2007a) used 
this multiplicity function to conclude that there are fewer 
LRG satellites of other LRG galaxies than predicted from 
N-body simulations. Furthermore, White et al. (2007) 
noted that the apparent lack of evolution in the cluster- 
ing strength of massive galaxies since 2^1 implies that 
these galaxies themselves must be merging as the under- 
lying dark matter distribution has undergone substantial 
merging during that epoch. Wake et al. (2008) measure 
the evolution of LRG clustering from z = 0.55 to z = 0.2 
and find that it is consistent with the idea that LRGs 
which originally resided in different halos merged to cre- 



ate a single galaxy when their host haloes merged. From 
the measured clustering of red galaxies in the NDWFS 
Bootes field. White et al. (2007) estimate that 1/3 of the 
LRGs which are satellites galaxies of another LRG have 
merged or been destroyed between z = 0.9 and z = 0.5. 

One model suggested to explain the deficit of LRG 
satellites suggests that the stars from late mergers onto 
massive galaxies feed the growth of an intracluster-light 
(ICL) type of extended envelope rather than the central 
galaxy. Conroy et al. (2007b) recently simulated the dis- 
sipationless evolution of galaxies since z = 1 and find 
that a model in which > 80% of the stars from merged 
satellites go into a low surface brightness extended stel- 
lar halo such as an ICL best predicts measurements of 
the galaxy stellar mass function and the observed distri- 
bution of ICL and brightest cluster galaxies in the local 
universe. If the total stellar content of the most mas- 
sive haloes grow considerably at z < 1 but the accreted 
stellar content resides in an extended, diffuse, envelope 
around the central galaxy, the total luminosity function 
of massive galaxies as measured by our technique would 
remain unchanged. 

It is clear from our observations that massive red galax- 
ies evolve in a systematically different manner than L* 
red galaxies. While the stellar mass in L* red galax- 
ies has doubled since z = 1, our analysis implies the 
mass in the L > 3L* red galaxies has grown, at most, 
by 50% over the same epoch. The growth of clusters 
and groups, including the intracluster light, may play a 
role in shaping the massive end of the red galaxy mass 
function while the lower-mass red galaxies are formed 
through the quenching of star forming galaxies at low 
redshifts. Alternatively, if the processes that govern star 
formation at the epoch of massive red galaxy formation 
are systematically different from those which govern star 
formation at z < 1, our analysis may underestimate the 
number density evolution in our sample. In the follow- 
ing section, we explore the impact that an evolving IMF 
would have on our analysis. 

4.4. Implication in the Presence of an Evolving Initial 
Mass Function 

Throughout all of our analyses, the slope of the stellar 
IMF is held fixed. While our dataset is not sufficient to 
constrain any evolution in the IMF of massive galaxies, 
if this evolution exists, it can strongly affect our con- 
clusions. Local measurements of the IMF show that at 
M > IMq the IMF follows a power-law {M/Mq cx M"^; 
X = 1.3) with a turnover at lower masses (Salpeter 1955; 
Kroupa 2001; Chabrier 2003). For this discussion, we will 
only consider the IMF at M > IMq; lower-mass stars, 
while contributing significant stellar mass to the galaxy, 
do not contribute significantly to the galaxy luminosity 
and thus play a negligible role in the evolution of the M/L 
ratio compared to variations in more massive stars. Sug- 
gestions of top-heavy IMFs have been found in environ- 
ments dominated by violent star-formation (Rieke et al. 
1993; McCrady et al. 2003; Figer et al. 1999; Stolte et al. 
2005; Maness et al. 2007). Also, one may expect the IMF 
to evolve with redshift as the temperature of the cosmic 
microwave background begins to dominate over tempera- 
tures typically found in Galactic prestellar cores (Larson 
1998). Recently, van Dokkum (2007) compared the lumi- 
nosity evolution of galaxies in clusters at 0.02 < z < 0.83, 
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coupled with the color-evolution of these systems, to test 
models of IMF evolution in early-type galaxies. These 
data prefer a logarithmic slope of x = — O.S^q ^, consid- 
erably flatter than x = 1.3 derived in the Milky Way 
disk. Similarly, Dave (2007) used hydro-dynamical mod- 
els of galaxy formation and observations of the correla- 
tion between galaxy stellar mass and star formation rate 
to = 2 to suggest that the characteristic mass at which 
the IMF turns over, M , evolves strongly with redshift : 

To explore the importance of the assumed IMF slope 
on the inferred density evolution in the LRG population, 
we show luminosity evolution tracks predicted using the 
fits of van Dokkum (2007) for SSPs formed at z = 2.0 
and z — 6.0 in Figure 16; the _B-band luminosity evolu- 
tion in each of the three tracks has been normalized to 
z = 0.3. The details of these models can be found in van 
Dokkum (2007). Briefly, these tracks show the expected 
luminosity evolution given three different IMF slopes us- 
ing Maraston (2005) synthesis models and [Fe/II]=:0.35. 
For slopes shallower than x = 1.3, our current passive 
evolution correction will systematically undercorrect for 
the passive fading of stars which will lead to significant 
underestimations of the density evolution experienced by 
these galaxies. For example, if we underestimate the lu- 
minosity evolution from z = 0.8 to z = 0.3 by 0.2 mag, 
we would conclude that the massive galaxy population 
has evolved little since z = 0.8 when, in actuality, the 
number density of these massive systems has grown by 
a factor of two. Clearly, more detailed constraints are 
needed on the fraction of high mass to low mass stars 
in these galaxies in order to place any evolutionary mea- 
surement into proper context. 
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Fig. 16. — B-band luminosity evolution based on initial mass 
functions with different slopes using the fits presented in van 
Dokkum (2007). The grey lines show the expected evolution of 
an SSP formed at z = 6 while the black lines show the trends for 
z = 2; all of the tracks have been normalized at z = 0.3. If galax- 
ies in our sample have IMF slopes shallower than the traditional 
X = 1.3 Salpeter (1955) value, we would underestimate the evolu- 
tion of galaxies at z; = 0.8 by > 0.15 mag by utilizing synthesis 
models based on the Salpeter (1955) IMF. 



4.5. Measurements of Massive Galaxy Luminosity 
Functions Using Aperture Luminosities 

Comparisons of several recent studies of the evolution 
of the red galaxy luminosity function since z ^ 1 have re- 
vealed a number of possible systematic differences which 
have been attributed to differences in the methods used 
to measure the total galaxy luminosities. For example. 
Brown et al. (2007) find that the stellar mass of the red 
galaxy population has grown by of a factor of 2 since 
z = 1.0 while results from DEEP2 suggest growth of a 
factor of 4 during the same epoch (Willmer et al. 2006; 
Faber et al. 2007). One alternative is to measure the 
luminosity of each galaxy in an aperture of fixed phys- 
ical size and to study the evolution of the luminosity 
function based on this quantity. This method removes 
the systematics introduced by comparing analyses done 
with fixed angular size aperture or extrapolations to the 
total galaxy flux. Furthermore, extrapolations to a total 
brightness requires careful treatment of the low surface- 
brightness outer isophotes which are quite difficult to 
photometer without very deep imaging. It is important 
to note, however, that the evolution of the luminosity 
within a fixed physical aperture size addresses a slightly 
different question than the total luminosity function; in- 
stead of tracking the total contribution of starlight, we 
instead focus on the growth of the stellar mass only in 
the inner region of the galaxy. Depending on the physical 
aperture size chosen, these luminosity measurements will 
not only be affected by the total starlight in the galaxy 
but also by the central concentration. Furthermore, the 
aperture luminosity function and total luminosity func- 
tion may exhibit different evolution if the ratio of the lu- 
minosity within the physical aperture to the total galaxy 
luminosity changes with time. For example, the aperture 
to total luminosity ratio may change if significant mass is 
accreted at large radii or the stellar concentration evolves 
due to recent merger activity. 

To investigate this method, we measure the evolution 
of the luminosity within the inner 2{)h~^ kpc for each 
galaxy in our sample. We choose 20/i~^ kpc radii aper- 
tures as this size will enclose a majority of the galaxy 
light, thus minimizing the effects on color gradients and 
galaxy concentration on our results, and yet not be too 
large such that the photometric errors due to sky sub- 
traction uncertainties become significant. For the low- 
redshift SDSS galaxy samples, we make use of the mea- 
sured aperture fiuxes at fixed angular sizes output by the 
SDSS pipeline. For reference, the SDSS pipeline mea- 
sures galaxy flux in apertures with radii of 0.23, 0.68, 
1.03, 1.76, '3.0, 4.63, 7.43, 11.42, 18.20, 28.20. 44.21, 
69.00, 107.81, 168.20, and 263.00 arcseconds (see Table 
7 in Stoughton et al. 2002). Based on the measured red- 
shift of each galaxy in our sample, we interpolate the 
measured aperture photometry to the radius correspond- 
ing to 20/i~^ kpc at the redshift of the galaxy. In order to 
measure the fluxes of our z ^ 0.9 galaxies at the highest 
possible signal-to-noise, we photometer these galaxies di- 
rectly from the SDSS imaging data. As our high-redshift 
sample was constructed from galaxies lying in the SDSS 
Southern Survey region, which has been scanned sev- 
eral times over the course of the survey, we construct a 
coadded image of 90/i~^ kpc x 90/i~^ kpc around each 
of our sample galaxies. Only data with seeing less than 
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Fig. 17. — Evolution of the luminosity function based upon 
luminosities contained within the central 20/i~^ kpc of massive 
galaxies. No significant differences are seen when the evolution of 
this central flux compared to the total galaxy luminosity functions 
presented in Figure 9. Measuring luminosities in apertures of fixed 
physical size eliminates systematic differences in estimates of the 
total galaxy fiux and thus will allow for more robust comparisons 
between future samples. 



1.5 arcseconds was used to construct the postage stamps. 
Before coadding each of the individual SDSS frames, we 
do not account for the seeing variations between each 
run; this has a negligible effect on the aperture fluxes on 
the scales we consider here. On each coadded postage 
stamps, known sources were masked out to avoid con- 
tamination and the flux of each galaxy was measured in 
a 20h^^ kpc radius aperture. 

Figure 17 and Tables 6 and 7 show the aperture magni- 
tude luminosity functions as a function of redshift. The 
aperture luminosity functions shown in Figure 17 show 
some systematic differences compared to the total lumi- 
nosity functions presented in Figure 9. At fixed lumi- 
nosity, the aperture luminosity function reports a sys- 
tematically smaller number density than the total lumi- 
nosity function. As the aperture luminosity measure- 
ments do not measure the full galaxy flux (with a me- 
dian Mapcr — Aftotai '^0.15 mag), the aperture luminos- 
ity function is shifted toward fainter magnitudes com- 
pared to the total luminosity function. Secondly, the 
number density falls off more rapidly toward more lumi- 
nous galaxies when aperture magnitudes are considered 
rather than total luminosities. This appears to be due 
to differential aperture losses as a function of luminosity; 
more luminous early-type galaxies have larger effective 
radii and thus more flux is missed by a fixed physical 
size aperture. While the shape and normalization of the 
aperture luminosity function have systematic differences 
with the total luminosity function, the aperture luminos- 
ity functions show little evolution in the 0.1 < z < 0.9 
range after the effects of passive evolution are removed 
just as is seen for the total galaxy luminosity function. 

The squares on Figure 11 and values in columns (3) 
and (5) of Table 4 and columns (3) and (6) of Table 5 
show the lack of evolution quantitatively - while the lumi- 



nosities computed using physically sized apertures were 
systematically fainter than the total galaxy luminosities, 
as expected, the evolution of the central 20 kpc of 
these massive red galaxies appears to follow the evolu- 
tion of the ensemble starlight. These measurements can 
provide a benchmark for future comparisons of the lumi- 
nosity function without the need to correct for systematic 
differences between the photometric methods used. 

5. SPECTRAL EVOLUTION OF MASSIVE GALAXIES 
SINCE Z ~ 0.9 

While each of our individual MMT galaxy spectra have 
too low signal-to-noise to perform any detailed measure- 
ments of line strengths, averaging the entire sample re- 
sults in a modest quality spectrum which can be used to 
measure the change in the spectral structure of massive 
red galaxies since z ^ 0.9. We construct the average 
LRG spectra in each redshift bin used to calculate our 
luminosity functions presented above : 0.1 < z < 0.2, 
0.2 < 2 < 0.3, 0.3 < z < 0.4, and 0.7 < z < 0.9. 
We limit the luminosity of the galaxies used in this 
analysis to the evolution-corrected magnitude range of 
—23 < Mo. 3^ — 51og h < —22 to focus on galaxies for 
which we are very complete. After masking within 10 
A of each of the strong emission lines arising from the 
Earth's atmosphere, we shift the observed spectrum of 
each galaxy to the rest-frame and normalize it by the av- 
erage fiux between 4100-4200A. We construct the mean 
spectrum by weighting each individual spectrum with the 
same weight assigned to that galaxy when calculating the 
luminosity function (including the l/Knax) and thus con- 
struct the composite spectrum of a typical galaxy in each 
of our redshift bins. 




3600 3800 4000 4200 4400 4600 
Wavlength (A) 

Fig. 18. — Average spectrum of LRGs since z = 0.9. Each 
composite spectrum shows features characteristic of old stellar 
populations while the highest redshift spectrum shows enhanced 
[O II]A3727 emission and stronger Balmer absorption indicating 
the presence of younger stars. The location of Balmer features are 
marked by vertical bars. As discussed in §5, we model the high- 
redshift average spectrum with a passively faded version of the 
low-redshift composite combined with a recent frosting of young 
stars. We find at most 5% of the stellar mass in the average high- 
redshift LRG has formed within IGyr o{ z = 0.9. 



16 



Cool et al. 



TABLE 6 

LHO 2()/i~^KPC Ai'EiiTrHE Lr.MiNDsrn' FrxcxiDXS with No Evou-jiox Cohhectiox 



logio Galaxy Number Density'' logio Galaxy Number Density'' 



1^-5 logh 


0.1<z< 0.2 


0.2 < 2 < 0.3 


0.3 < z < 0.4 


Mo.3^ - 5 log ^ 


0.7 < z < 0.9 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


-21.55 


-3.75 ±0.01 


-3.72 ±0.01 


-3.49 ± 0.04 


-21.59 


-3.27 ±0.07 


-21.65 


-3.92 ±0.01 


-3.84 ±0.01 


-3.59 ±0.03 


-21.77 


-3.33 ± 0.07 


-21.75 


-4.09 ±0.02 


-3.99 ± 0.01 


-3.73 ±0.02 


-21.95 


-3.51 ± 0.08 


-21.85 


-4.30 ±0.02 


-4.12 ± 0.04 


-3.88 ±0.04 


-22.13 


-3.74 ±0.09 


-21.95 


-4.55 ±0.03 


-4.39 ± 0.01 


-4.08 ±0.05 


-22.31 


-4.01 ±0.10 


-22.05 


-4.93 ±0.05 


-4.62 ± 0.02 


-4.33 ±0.01 


-22.49 


-4.53 ± 0.15 


-22.15 


-5.15 ±0.06 


-4.94 ± 0.03 


-4.58 ±0.02 






-22.25 


-5.64 ±0.12 


-5.23 ± 0.04 


-4.85 ±0.02 






-22.35 


-5.93 ±0.16 


-5.52 ± 0.05 


-5.18 ±0.02 






-22.45 


-6.30 ±0.25 


-6.19 ± 0.12 


-5.53 ±0.03 






-22.55 


-6.45 ± 0.43 


-6.63 ± 0.19 


-5.86 ±0.04 






-22.65 


-6.64 ±0.43 


-6.63 ± 0.19 


-6.20 ±0.06 







''All number densities are expressed in units of Mpc ^ Mag 



TABLE 7 

LRCt 20/i~^ki'<' Ai'i;irrrnE LrMixosriT Erxciioxs Aftek I'assixx livoLrxiox C'oKiiEtTiox 



logio Galaxy Number Density" logio Galaxy Number Density" 



— 5 logh 


0.1<z< 0.2 


0.2 < 2 < 0.3 


0.3 < 2 < 0.4 


Mo. 3^ - 5 log ^ 


0.7 < 2 < 0.9 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


-21.55 


-3.51 ±0.01 


-3.44 ± 0.02 


-3.51 ±0.02 


-21.59 


-3.61 ±0.08 


-21.65 


-3.62 ±0.01 


-3.56 ± 0.02 


-3.62 ±0.02 


-21.77 


-3.84 ± 0.09 


-21.75 


-3.77 ±0.01 


-3.74 ±0.02 


-3.77 ±0.02 


-21.95 


-4.12 ± 0.11 


-21.85 


-3.91 ±0.02 


-3.93 ± 0.02 


-3.93 ±0.05 


-22.13 


-4.51 ± 0.15 


-21.95 


-4.13 ±0.02 


-4.16 ± 0.01 


-4.19 ±0.01 


-22.31 


-5.09 ± 0.22 


-22.05 


-4.33 ±0.02 


-4.37 ±0.02 


-4.42 ±0.01 


-22.49 


-5.81 ±0.23 


-22.15 


-4.64 ±0.03 


-4.62 ± 0.02 


-4.69 ±0.02 






-22.25 


-4.92 ±0.05 


-4.98 ± 0.03 


-4.98 ±0.02 






-22.35 


-5.26 ±0.07 


-5.27 ±0.05 


-5.29 ±0.03 






-22.45 


-5.65 ±0.12 


-5.69 ± 0.07 


-5.72 ±0.06 






-22.55 


-6.31 ±0.25 


-6.31 ± 0.15 


-6.19 ±0.09 






-22.65 


-6.72 ±0.43 


-6.62 ± 0.22 


-6.95 ±0.21 







"All number densities are expressed in units of h'^ Mpc ^ Mag 



Figure 18 shows the coaddcd spectra of massive red 
galaxies from z = 0.1 to z = 0.9. Each of the compos- 
ite spectra look quite similar showing the strong spectral 
features characteristic to old stellar populations. While 
the high-redshift composite spectrum clearly shows en- 
hanced [O II] emission compared to the lower-redshift 
spectra other differences between the spectra are more 
subtle. Figure 19 shows the measured H5 and G-band 
at 4300 A absorption equivalent width, from our compos- 
ite spectra. A solar-mctallicity stellar population formed 
at z = 2 using a Salpeter (1955) IMF with subsequent 
passive fading is shown with the solid line. Our mea- 
surements are broadly consistent with the passive fading 
of stars since z ~ 0.9. Note that we make no claim that 
since these points lie near the solar-metallicity track that 
we expect these galaxies to have solar metallicity or have 
a given age. It has been shown (e.g. Eisenstein et al. 
2003; Cool et al. 2006) that LRGs show a-enhancements 
compared to solar and also that the age and metallic- 
ity of the stellar populations one might derive from most 
spectral indicies are degenerate. Instead, we simply il- 
lustrate that the data follow the same trend expected for 
a passively fading population. 

In order to model the amount of recent star forma- 



tion activity allowed by our high-redshift composite spec- 
trum, we model it as the linear combination of a passively 
faded version of our low-redshift spectrum plus a frost- 
ing of more recent star formation activity. The lowest- 
redshift composite is well fit by a 7.0 Gyr, solar metallic- 
ity, population. Thus, we model our high-redshift com- 
posite as the non- negative linear sum of a 1.9 Gyr popula- 
tion - the universe has aged by 5.1 Gyr between z = 0.8 
to z = 0.15 - and a frosting of either lOMyr, lOOMyr, 
or IGyr stars. We find that the high-redshift composite 
is best modeled by a single-age population at 1.9 Gyr 
with no need for the presence of younger stars save for 
the [O ll] which may be generated by either young stars 
or enhanced AGN activity. We can constrain the pres- 
ence of lOMyr, lOOMyr, and IGyr stars to contribute less 
than 0.1%, 0.5%, and 5% of the stellar mass based on our 
spectral fits with 99% confidence. Thus, it appears that 
high-redshift LRGs have enhanced signatures of youth 
compared to their low-redshift counterparts due to the 
passive evolution of their stellar populations. We find no 
signatures of more recent star formation activity in our 
high-redshift sample indicative of recent gas-rich mergers 
at z ~ 0.9. 

The evolution of the average spectrum presented here 
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may be underestimated in the event that galaxies with 
weak absorption lines are preferentially removed from 
the sample due to redshift determination failures. We 
do not expect our spectroscopy to be biased in this way, 
however. Primarily, as the absorption line strength is 
correlated with the total galaxy luminosity, we expect 
the galaxies with weak lines to have luminosities fainter 
than the limits imposed in creating our composite spec- 
tra. Secondly, we would expect the presence on [O II] 
emission to allow redshift determination even if the ab- 
sorption lines were very weak. To examine this effect, 
we refit each of our galaxies after masking out the wave- 
lengths affected by the [O II] emission line and find that 
only 3 of the galaxies in our sample had sufficiently weak 
absorption lines that the presence of [O ll] dominated the 
redshift fitting. 

6. CONCLUSIONS 

Massive galaxies serve as probes of the merger history 
of the universe as these galaxies have participated most 
heavily in the merger process. Using samples of mas- 
sive {L > 3L*) red galaxies observed by SDSS at low 
redshift augmented with a new spectroscopic sample of 
galaxies targeted from deep SDSS coadded photometry 
and observed with the MMT, we have measured the evo- 
lution of massive red galaxies at 0.1 < z < 0.9. Our 
sample is currently the largest collection of massive red 
galaxies spectroscopically observed at z ~ 0.9 and thus 
provides an excellent tool for constraining the evolution 
of the most massive galactic systems over half of cosmic 
history. 

After correcting for passive evolution using a non- 
evolving Salpeter (1955) IMF, we find the magnitude at 
which the integrated number density of the LRG pop- 
ulation has reached 10~^'^/i^ Mpc~^ is consistent with 
constant with a best-fit evolution of 0.03 ±0.08 mag from 





o 


0.1<z<0.2 ■ 






o.2<z<o.3 : 




□ 


0.3<z<0.4 ■ 




^\ O 


0.7<z<0.9 : 






A * 




1 . . . . 


1 . . \ . 



4.0 4.5 5.0 5.5 

G4300 (EW) 

Fig. 19. — Equivalent widths of US and G-band absorption 
features from the composite galaxy spectra. The style of data 
point corresponds to the redshift of the composite spectrum; star: 
0.1 < 2 < 0.2, asterisk: 0.2 < z < 0.3, square: 0.3 < z < 0.4, and 
circle: 0.7 < z < 0.9. The solid line shows the expected trend for 
a solar-metallicity galaxy formed at 2; = 2 from Bruzual & Char- 
lot (2003) models. Errors are comparable to the size of each data 
point. While only illustrative, the observed composite spectra show 
similar trends as that expected of a passively fading population. 



z — I to z = 0. Simple toy models for the merger his- 
tories of massive red galaxies indicate that 1:1 merger 
rates larger than 25% are disfavored at 50% confidence 
and merger rates larger than 40% are ruled out at 99% 
significance. Even if lower-mass mergers are considered, 
we find that the total stellar mass contained in massive 
red galaxies must not have grown by more than ^ 
since z = 0.9. This growth rate starkly contrasts the 
factor of 2-4 in stellar mass growth observed in L* red 
galaxies over the same epoch. The processes that regu- 
late the growth of massive red galaxies and yet allow the 
large growth observed in the L* red galaxy population 
are poorly understood. As the most massive galaxies re- 
side in group and cluster sized haloes, the processes that 
govern the assembly of clusters or the growth or an in- 
tracluster stellar envelope may play an important role in 
the shaping of LRGs. 

The evolution in the average LRG spectrum to high 
redshift also supports a purely passive fading of LRGs 
since z ~ 0.9. The composite spectrum of our high- 
redshift LRGs is well-described by a passively faded ver- 
sion of the average galaxy spectrum at0.1<z<0.2. No 
recent star formation is needed to explain our composite 
spectrum at z = 0.9; we constrain the mass fraction of 
lOMyr, lOOMyr, and IGyr stars to be less than 0.1%, 
0.5%, and 5% with 99% confidence. Star formation in 
these LRGs must have completely ended hy z ^ 0.9 and 
very few blue stars must have been accreted since that 
epoch. 

While our sample comprises the largest spectroscopic 
sample of massive red galaxies at z ^ 0.9 collected to 
date, a sample of 300 galaxies suffers from small num- 
bers of objects per luminosity bin, especially at the high- 
est masses. Future surveys aiming to collect spectro- 
scopic samples of many thousand LRGs at redshifts up 
to z ^ 0.7, while at slightly lower redshifts, will have 
the statistics to place tighter constraints on the overall 
density evolution of the massive red galaxy population 
as well as to study the evolution in the LRG luminosity 
function shape to constrain the role of mass-dependent 
processes which regulate LRG growth. 
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